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^ : Abstract 

, ^ , The final state of the decay — > vr^vr vr'^ is analyzed in terms of isospin eigenstates. It is 
shown that the final state is dominated by the isospin-0 component. This suggests that isospin 

■ considerations may provide insight into this and perhaps other D'^-meson decay. We also discuss 

r~| . the isospin nature of the nonresonant contribution in the decay, which can be further understood 

' \ by studying the decay tt^tt^tt^. 
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I. INTRODUCTION 



An analysis of the resonant sub-structure in the decay tt^tt tt^ was recently per- 



formed by the BABAR collaboration 



1[. The Dalitz-plot distribution of the D 



n n TT 



events (Fig. [T]) shows a clear six-fold symmetry, with the probability density function van- 
ishing along three axes. As first described by Zemach and noted in Ref. [jj, this behavior 
is indicative of a final state with isospin / = 0. 
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FIG. 1: The Dalitz-plot distribution of it~^tt tt^ events, from Ref. 



In the BABAR analysis, the Dalitz-plot distribution is described by a probability density 
function formed from a wave function taken to be the sum of contributions, 

Nr 

i^is+,S-) =^Brgr{s+,S^), (l) 
r 

where s+ = {pt,+ +Pno)'^ and s_ = {pt,- +j97ro)^ are the squared invariant masses of the 7r+7r° 
and 7r~7r° pairs, respectively, Br is a complex coefficient, and gr{s+,s^) is the distribution 
of contribution r, whose functional form is outlined in Ref. The definitions of gr{s+, s_) 
used here differ from that of Ref. [l| , in that we define these functions to be normalized over 
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the Dalitz plot, 

j ds+ds-\gr{s+,s^)f = 1. (2) 
The values for the Br coefficients consistent with Eqs. and ([2]) are reproduced in Table [H 

TABLE I: Amplitude coefficients = \Br \ e*'^'' of the contributing final states of the decay —i- 
vr-vr+vrO, adapted from Ref. Q. The /o(400) was labeled a(400) in Ref. 



Final state r 


Amplitude \Br\ 


Phase (pr (°) 


Nonresonant 


0.106 ± 0.013 ± 0.014 


-11±4±2 


p{770)+TT- 


1 


0.0 




0.588 ± 0.006 ± 0.002 


16.2±0.6±0.4 


p(770)-7r+ 


0.714 ± 0.008 ± 0.002 


-2.0±0.6±0.6 


p(1450)+7r- 


0.040 ± 0.011 ± 0.024 


-146±18±24 


p(1450)07r0 


0.062 ± 0.012 ± 0.007 


10±8±13 


p(1450)-7r+ 


0.154 ± 0.010 ± 0.007 


16±3±3 


p(1700)+7r- 


0.236 ± 0.019 ± 0.014 


-17±2±3 


p(1700)0^° 


0.267 ± 0.016 ± 0.014 


-17±2±2 


p(1700)-7r+ 


0.210 ± 0.012 ± 0.007 


-50±3±3 


/o(980)^o 


0.056 ± 0.005 ± 0.006 


-59±5±4 


/o(1370)7rO 


0.072 ± 0.010 ± 0.010 


156±9±6 


/o(1500)7rO 


0.074 ± 0.007 ± 0.007 


12±9±4 


/o(1710)7rO 


0.072 ± 0.010 ± 0.011 


51±8±7 


/2(1270)7r0 


0.130 ± 0.005 ± 0.026 


-171±3±4 


/o(400)^o 


0.104 ± 0.008 ± 0.017 


8±4±8 



The goal of this paper is to quantify the extent to which the 1 = component dominates 
the final state and learn about the contributions of the other isospin eigenstates. In Sec- 
tioniniwe perform an isospin analysis of the tc^tt^ti^ final state. The observed dominance of 
the 1 = component suggests that isospin considerations are more useful for developing an 
understanding of this decay. In Section [XTTl we discuss our results, the nature of the nonres- 
onant contribution to the decay, a possible mechanism for the observed / = dominance, 
and further measurements that will help clarify outstanding questions. 
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II. ISOSPIN DECOMPOSITION 



Next, we analyze the decay D'^ 7r"'"7r~7r'^ in terms of isospin eigenstates. The 3-pion final 
state can be described in terms of the total isospin /, the isospin I12 of two of the three pions, 
and the 2;-projection P, which is always for this final state. The seven eigenstates |/(/i2)) 
of these quantum numbers that also satisfy P = can be written as a linear combination 
of the three-pion final states using the appropriate Clebsh-Gordan coefficients: 



|3(2) 

|2(2) 

|1(2) 

|2(1) 

|1(1) 
|0(1) 

|1(0) 



^(1+0-) + |0 + -) + 1+ - 0) + I- + 0) + |0 - +) + 1-0+) + 2 |000)) 
^(|+0-) + |0 + -)-|0-+)-|-0+)), 

(I+O-) + |0 + -) + |0 - +) + 1-0+)) - 2(^1+ - 0) + I- + 0)) - 4 |000) 
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+0-) - |o + -) + 2(^1+ - 0) - I- + 0)) + |0 - +) - 1-0+) 

^(|+0-)-|0 + -)-|0-+) + |-0+)), 

(I+O-) - |0 + -) - 1+ - 0) + I- + 0) + |0 - +) - 1-0+)) , 

-^(|+-0)-|000) + |- + 0)), (3) 

where we have used the notation 

I+0-) = |l,l)|l,0)|l,-l) = |7r+)|7r°)|7r-), 
1000) = |l,0)|l,0)|l,0) = |7r°>|7r0>|7r°>, (4) 

etc., and it is implied that the first two pions are in an isospin eigenstate whose eigenvalue 
is indicated by the bracketed number I^. 

The three states in Eq. ([3]) for which J12 = 1 are identified as those with a p(770), 
p(1450), or p(1700). We denote these states as p„7r according to their radial excitation 
quantum number G {1, 2, 3}, and use p^, p^, and p^ to indicate any linear combination of 
these states with specific electric charge. We define the p states to be 



P 



l,l) = i=(|+0)-|0+) 



|P°> = -|i,o) = i=(|-+)-|+-) 

|P"> = |l,-l) = ^(|0-)-|-0)), (5) 



where the minus sign in the Ip*^) definition imphes that there is no sign change under cychc 
permutations of the three pions, maintaining consistency with the definitions used in Ref . [l| . 
Given Eq. (E]), the lu = 1 states in Eq. can be written as 

|2(1)) = -^(|p%->-2|pV> + |p--+>), 

|0(1)) = ^(|p+vr-> + |pV> + |p-7r+>), (6) 

where the sign of each Ipvr) state is such that it is symmetric under cychc permutations of 
the three pions and anti-symmetric under the exchange of any pair of pions. 

The 7r+7r^7r° part of the state |1(0)) is identified as the sum of the contributions involving 
the two-body, / = resonances fi, with i = 0,2. We therefore write 

|l(0)) = i=(V2|K>-|000)). (7) 

Since there are no / = 2 resonances in Table [H the /12 = 2 states in Eq. ([3]) have no res- 
onant contributions. However, the symmetry of the tt+tt^tt'^ components of |3(2)) indicates 
that it may be identified with the nonresonant contribution of Table HI Alternatively, it may 
constitute the 7T~^n^7T^ component of the symmetric 1 = 1 state 

= ^|1(2)) + ^|1(0)) 

^ f I+0-) + |0 + -) + 1+ - 0) + I- + 0) + |0 - +) + 1-0+) - 3 |000)) . (8) 
15V / 



as 

In principle, the observed nonresonant state may be a superposition of m^*)) and |3(2)). 
However, the |1(5')) state is expected to dominate, due to the following argument. The 
four-quark final state produced by the weak decay c — > ddu, shown in Fig. [2], cannot have 
/ = 3. Since production of the third qq pair will be dominated by the strong-interaction, 
it will not change the total isospin. Therefore, J = 3 is disfavored. It is also possible that 
a very broad, tt+tt^ S-wave resonance is present in these decays, and that it was partly 
described by the constant nonresonant term in the fit in Ref. In that case, it would 
contribute only to the |1(0)) isospin eigenstate. 

In what follows, we take the nonresonant contribution |A^-R) to be due only to I^S*)). 
Then Eqs. ([7]) and ([8]) yield the relation 



We now reorder the terms of Eq. ([T]) according to their I12 eigenvalues: 

V'(s+,s_) = 5nr5'nr(s+,s_) 

+ BpOj^o 5(p0^o(s+, S-) 

+ Bf^o gf^o{s+,s^), (10) 

where the first term is the nonresonant term, the last is a sum over the six final states with 
/12 = resonances listed at the bottom of Table [U and each of the second, third, and fourth 
terms is a sum over the three I12 = 1 pvr states. For example, 

S + 

gp+n-{s+,s_) = J exp [-i5p+^-] , (11) 

where 



Sp+n- = J2BpU-9pU-{s+,S-), 

n=l 

Sp+n- = arg(5'p+^-) , 



iVp+TT- = yy ds+ds_\Sp+^-\ , (12) 

and pn {n = 1,2,3) indicates the three p resonances of Table [H With these definitions, 
the wave function ^fp+Tr- (s+, s_) is explicitly normalized and has vanishing average phase. 
Requiring that Eq. ffTUl) be identical to ([T]) leads to the following values for the coefficients 
of Eq. (Uni): 





= 0.1066 e-*^^-^°. 




Bp+TT- 


= Np+^- exp [i5p+^~ 


] = 1.1976 e-'^-=^° 


BpO^O 


= NpO^^o exp [i^pOjro] 


= 0.8867 e^^-=^°. 


Bp-n+ 


= Np-^+ exp [i6p-n+ 


] = 1.0077 e-'^-2° 


BfjrO 


= NfT^o exp [«5/7ro] = 


: 0.0700 e*^°°°. 



(13) 

where the symbols Ns and 6s for final state s are defined analogously to Eq. f|T2l) . The value 
of Snr is taken from Tableland the rest are calculated numerically as in Eqs. (ITT]) and (fT2l) . 
The phase convention is that of Table [H namely, 6^+^- = 0. 



Next, we write the wave function of Eq. (fTOj) as a sum over the Dahtz-plot representations 
of the eigenstates of I and /12 of Eq. ([3]): 



= Ci(2)Mi(2)(s+,S_) 
+ C2(1)M2(1)(S + ,S„) 

+ Ci(i)Mi(i)(s+,s_) 

+ ^0(1)^-4(1) (-5+, S_) 

+ Ci(o)Mi(o)(s+,s_), (14) 

where M7(/^2)('5+) ■5-) is the normahzed distribution function of the eigenstate |/(/i2)), ob- 
tained by hnearly combining the functions gx{s+., sJ) of Eq. ffTOj) with the coefficients of 
either Eq. ([6]), ([7]), or ([9]). Terms for |3(2)) and |2(2)) were not included in Eq. (IT^ . as rea- 
soned earher. Then from the definition of M/(/j2)(s+, s_) follows the desired transformation 
between the resonance-based fit coefficients and the isospin coefficients: 



Cl(2) = 


2 -DNR, 




C2{1) = 




— 2Bp0^o + -Bp-TT 


Cl(l) = 


7! 




Co(i) = 


7! 


+ -BpO^O + -Bp-7r+ 



'^i(o) = \I^Bf^o + ^^C^2), (15) 

where the expressions for Ci(o) and Ci(2) were chosen so as to satisfy the n^n^n^ projection 
of Eqs. (ED and (ED. 

Taking the numerical values of the Br coefficients from Eq. ( |T3D and Table [H Eq. ( ITSP 
gives 

Ci(2) = (0.0629 ± 0.0028) exp [i (-8.9 ± 2.6)°] , 
C2(i) = (0.1395 ± 0.0016) exp [i (-42.5 ± 0.7)°] , 
Ci(i) = (0.0814 ±0.0023) exp (18.0 ±2.0)°], 

Co(i) = 1, 

Ci(o) = (0.0954 ± 0.0052) exp [i (14.5 ± 2.4)°] , (16) 
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where we have normahzed the coefficients so that Co(i) = 1- The errors reflect the full error 
matrix of the results presented in Tabled ^. The correlation matrix for these coefficients 
are given in Table [TTl 

Eq. flTB]) quantifies the observation, made qualitatively in Ref. |l| on the basis of the 
symmetry exhibited by the Dalitz-plot distribution, that the final state of the decay — > 
7r+7r~7r° is dominated by an / = component. 



TABLE II: Correlation matrix for the C/(7i2) amplitude coefficients of Eq. p!6l 





|Cl(2)l 


arg(Ci(2))| 


1^2(1)1 


arg(C2(i)) 


|Ci(i)l 


arg(Ci(i)) 


|Ci(o)l 


arg(Ci(o)) 


|Cl(2)l 


1 


-0.120 


0.105 


-0.018 


0.631 


0.110 


0.279 


0.657 


arg(Ci(2)) 


-0.120 


1 


0.062 


0.106 


-0.211 


0.539 


-0.760 


0.136 


1^2(1)1 


0.105 


0.062 


1 


0.008 


0.179 


0.029 


-0.017 


0.078 


arg(C2(i)) 


-0.018 


0.106 


0.008 


1 


0.148 


0.333 


0.110 


0.151 


|Ci{i)l 


0.631 


-0.211 


0.179 


0.148 


1 


0.050 


0.259 


0.288 


arg(Ci(i)) 


0.110 


0.539 


0.029 


0.333 


0.050 


1 


-0.296 


0.097 


|Ci(o)l 


0.279 


-0.760 


-0.017 


0.110 


0.259 


-0.296 


1 


0.077 


ai-g(Ci(o)) 


0.657 


0.136 


0.078 


0.151 


0.288 


0.097 


0.077 


1 



III. DISCUSSION AND CONCLUSIONS 

We have analyzed the relative contributions of different components to the decay — »• 
7r+7r^7r° using results published by BABAR [1]. It appears that isospin considerations may 
form a solid basis for understanding the observed decay pattern, as the amplitude of the 
1 0(1)) final state dominates by factors of seven or more over the other isospin components. 
This dominance has no natural explanation in the decay mechanisms suggested by the 
factorization-motivated diagrams of this decay, shown in Fig. [21 While factorization is 
useful in predicting the behavior of i?-meson decays, it is not as successful when applied 
to the lighter D mesons. The observed |0(1)) dominance in the decay tt+tt^tt'^ may 

lead to a better general understanding of charmed meson decays. Alternatively, perhaps the 
1 = component is enhanced by the presence of a yet-unknown and possibly broad state 
with this quantum number, which couples strongly to three pions. An inclusive search for 
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FIG. 2: Feynman diagrams for the decay vr^vr vr". With curly brackets indicating a res- 

onance, the diagrams correspond to the decays (a) — > p^tt~, (b) — > tt^/?^, and (c,d) 



such a state may answer this question. 

In conducting the isospin analysis, we took only the tt+tt^tt^ projections of the isospin- 
eigenstates |1(2)) and |1(0)). The CLEO collaboration 5| has set an upper limit of 3.4x 10~^ 
on the branching fraction B{D^ TT^Tr^iv^). Together with the BABAR [g] measurement of 



vr TT TT 



0^ 



1^1.493 ± 0.057)%, this implies an upper limit on the amplitude ratio 



A{D^ ^ Ti^TT^TT^) /A{D^ ^ 71^71^71^) < 0.15, cousistcut with the suppression seen in the 
coefficients Ci(2) and Ci{o), and the expectation from Eqs. ([7]) and 

As discussed above, the tt+tt^tt^ nonresonant amplitude may be a combination of |3(2)), 
mS*)), and a broad vr+vr" resonance term in |1(0)). If it is due only to the |3(2)), Eq. ([3]) pre- 
dicts the ratio between the nonresonant vr'^Tr'^Tr*^ and n^n^n^ amplitudes to be Rnr = a/2/3. 
By contrast, |l(S'))-dominance leads to Rnr = a/3/2, from Eq. ([8]). In the |1(0)) case, the 
ratio between the nonresonant vr^Tr^Tr*^ amplitude and the sum of the Jti^ and nonresonant 
7r+7r^7r° amplitudes should be 1/v^. We note that the ratio Rnb. = Vl-556 ± 0.012 is ob- 
served in Kl decays to three pions, where the nonresonant contribution accounts for over 
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95% of the branching fractions. The same situation exists in the decay 77 — > n'^n~n'^. This 
strengthens the justification of our choice to identify the nonresonant contribution with the 
1 1(5')) state. In any case, the arguments given here demonstrate that a measurement of the 
branching fraction B{D^ Tr^Tr^Tr") and, possibly, an analysis of this mode's Dalitz-plot dis- 
tribution should shed more light on the role of isospin symmetry in decays to three-pion 
final states. 
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